Geometrical enhancement of the electric field: 
Application of fractional calculus in nanoplasmonics 
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We developed an analytical approach, for a wave propagation in metal-dielectric nanostructures 
in the quasi-static limit. This consideration establishes a link between fractional geometry of the 
nanostructure and fractional integro-differentiation. The method is based on fractional calculus and 
permits to obtain analytical expressions for the electric field enhancement. 
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I. INTRODUCTION 

An interplay between the light and fractional nanos- 
tructures leads to a local giant enhancement of the elec- 
tric field. This is one of the main effect of nanoplasmon- 
ics [l|, Ql that concerns with the electromagnetic field on 
a nanoscale less then the wavelength inside a composite 
matter. There is a vast literature starting from seminal 
results on optics in fractal smoke [3j with further develop- 
ing, where strong fluctuations of a local electromagnetic 
field are taken into account, see detailed discussion in 
[J| and recent reviews [E S Q- ^ suitable mechanism 
of these phenomena have been proposed in many stud- 
ies [l|, y-Q)' s-'^d it relates to the optical properties of 
metal nanoparticles, and, mainly, is due to their abil- 
ity to produce giant and highly localized electromagnetic 
fields inside composite nanostructure. We are speaking 
about a localized surface plasmon resonance (SPR) (see, 
for example Ref. [i|, Q). Localized surface plasmons are 
charge density oscillations confined to conducting (metal- 
lic) nanostructures. A strong enhancement of the electro- 
magnetic field is possible if its frequency is near the fre- 
quency of these oscillations. An analogous phenomenon, 
also known in optics, is a so-called enhanced Raman scat- 
tering. This remarkable phenomenon has been firstly ob- 
served more than thirty years ago, see e.5., [a, 01 and 
properties of the Raman spectroscopy, has been used for 
a variety of applications. Theoretical explanation fol- 
low these experiments VA and still attract much attention 

In this paper we suggest an analytical consideration, 
based on fractional calculus, purposely to obtain an an- 
alytical expression for the electric field enhancement in 
nanoplasmonics. This approach establishes a relation be- 
tween fractal geometry of the nanostructure and frac- 
tional integro-differentiation. It is based on averaging an 
extensive physical value expressed by means of a smooth 
function over a fractal set that leads to fractional in- 
tegration fs]. We suggest a coarse graining procedure 
of the electric field in the Maxwell equation to treat a 
charge density term, which diverges everywhere in the 
case of fractional metal-dielectric composite media. This 
smoothing procedure makes possibility to obtain an equa- 
tion for an electric field in a closed form. Our approach 
is based on a seminal result [9|, where this scheme has 



been suggested for a Cantor set. In its eventual form, it 
has been presented in Ref. [10| (in Ch. 5). The main 
idea of embedding or filtering matter inside a fractal is 
the construction of a convolution integral 



M{x) = W{x)i. f{x) 



W{x-v)f{y)dy, (1) 



where the function W{x) obeys the following scaling re- 
lation W{x) = \W{bx) with a solution W{x) = i"A(x), 
where a — ^ and A{x) — A{hx) is a log-periodic func- 
tion with a period In h. When this scaling corresponds to 
a Cantor set with a — 2, and A(x) is defined explicitly, 
one performs averaging over this period [10| and obtains 
the convolution integral in the form of a fractional inte- 
gral 



{M{x)) = 



r(a) 



[x ~ x')°'-^f{x')dx' , 



(2) 



where V{a) 



-I + q/2 

Till 



and T{a) is a gamma function. 



One should recognize that the coarse graining procedure 
due to the averaging over the period In b is important for 
the application of fractional calculus, and "a bond be- 
tween fractal geometry and fractional calculus" can be 
established on some coarse graining geometry [l3|. This 
mathematical construction is relevant to study of electro- 
statics of real composite structures in the framework of 
coarse-graining Maxwell's equation. The main objective 
of the preset research is an analytical derivation of the 
electric field in a composite dielectric, which is subject to 
an external high-frequency electric field at the conditions 
when the quasistatic limit is valid. 



II. FRACTIONAL CALCULUS BRIEFLY 

Extended reviews of fractional calculus can be found 
e.g., in |ll| - [l3| . Fractional integration of the order of a 
is defined by the operator 



j:i(x) 



1 



r(a) 



f{y){x-yr-Hy, 



where a > 0, x > a. Fractional derivation was developed 
as a generalization of integer order derivatives and is de- 
fined as the inverse operation to the fractional integral. 



Therefore, the fractional derivative is defined as the in- 
verse operator to a-f", namely aD"f{x) = alx" f{^) ^^^ 
al" — aD~°'. Its exphcit form is 



aD^fix) 



1 



Ti-a) 



fiy)ix~yr'-'^dy. 



For arbitrary a > this integral diverges, and, as a result 
of this, a regularization procedure is introduced with two 
alternative definitions of aD". For an integer n defined 
asn — 1 < a < n, one obtains the Riemann-Liouville 
fractional derivative of the form 



^D^fix) = 4"-/(x) 



dx^ 



and fractional derivative in the Caputo form 



<;^D^fix)=j:-'^—f{x) 



dx^ 



(3) 



(4) 



There is no constraint on the lower limit a. For example, 
when a = 0, one has ^^D'^x^^ = ""^ff+i^^^^ , a, /3 > 0. 
This fractional derivation with the fixed low limit is also 
called the left fractional derivative. Another important 
property is D"/^ = /^~", where other indexes are omit- 
ted for brevity's sake. A convolution rule for the Laplace 
transform for < a < 1 

C[i:f{x)] = s-V~{s) (5) 

is commonly used in fractional calculus as well. 



III. QUASI-STATIC ELECTROMAGNETIC 
FIELD 



We consider a sample of a size R which consists of a 
fractal metal nanosystem embedded in a dielectric host 
medium. The system is subject to an external electro- 
magnetic field Eo(t) at an optical frequency uj, and we 
assume that the fractal-host inhomogeneities are much 
smaller than the light wavelength A = 2ttc/uj: Iq <C A, 
where c is the light speed and Iq is a minimal size of 
the self-similarity of the fractal volume. In this case a 
quasi-static approximation is valid |14| . Therefore, the 
nanosystem contains a metallic fractal with a fractal vol- 
ume Vd ~ R^ and the permittivity em(aj), which de- 
pends on optical frequency w, while the dielectric host 
has the volume Vh and the permittivity Sd- Following 
Refs. [a, LZ|, the wave magnetic component is not consid- 
ered, since the magnetic permeability /i = 1 in both me- 
dia. The electric component satisfies the Maxwell equa- 
tion in the Fourier frequency domain 



V- e(r,w)E(r,w) =0 



(6) 



The permittivity can be expressed by means of a char- 
acteristic function x(r) (a-Q in the form e(r, w) = 
£m(w)x(r)-|-ed[l— x(r)], where the characteristic function 



inside the fractal is x(r) = 1, r{x, y, z) €Vd, while inside 
the dielectric host it reads x(r) = 0, r{x,y,z) G Vh- At 
the boundaries the electric field consists of the external 
field and an induced field due to polarization of the frac- 
tal nanostructure. Splitting the electric field into the two 
components E(r) — E(r)-|-Ei(r), we looking for the reac- 
tion of the nanosystem on the external field, where E(r) 
is the electric field induced by Eq in the homogeneous 
nanosystem, when x(r) = for V r, or x(r) = I for V r 
and VE(r) = 0. Here Ei(r) results from inhomogeneity 
of the nanostructute. Thus Eq. ^ can be rewritten in 
the form 

x(r)V-Ei(r) - <zHV • Ei(r) + Ei(r) • Vx(r) 

= -E(r).Vx(r), (7) 



is a spectral parameter I5| . The 



where q(uj) = ^^-^-^ 

boundary condition for the normal component of the di- 
electric shift is D„(r G 5) = E„(r G 5), where S means 
the sample boundaries. 

The main complication of treating Eq. ([7]) is the po- 
larization charge density term ^ Vx(r). Since the char- 
acteristic function x(r) is discontinuous, Vx(r) diverges 
everywhere. Therefore, the next step of our considera- 
tion is coarse-graining the electric field, which is averag- 
ing Maxwell's Eq. ([T]). This procedure relates to inte- 
gration with the characteristic function J x(i")^Ei(r)(iF 
and evaluation of the integral / Vx(r) • E(r)dV^. To that 
end, let us consider a spherical volume of the radius 
r, such that /q ^ r ^ A. At this scale the magnetic 
component is not important and cannot be considered. 
In the sequel we will work with dimensionless variable 
r//o — > r. The electric field does not change at this scal- 
ing. A fractal mass inside the volume is M.{r) ~ r^, 
where < D < 3. Therefore, an average density of the 
metallic phase is of the order of r^~^. For the filtering 
mass inside the fractal, we consider convolution of Eq. 
dl]). We also, reasonably, suppose that the random frac- 
tal composite is isotropic, i.e., the (averaged) similarity 
exponents coincides along all directions. This yields that 
one takes into account only radius r and any changes 
in the inclination and azimuth angle directions can be 
neglected. Then, we have for the divergence 



V • Ei(r) = VrEr,i{r) = V^^i(r) . 



(8) 



Now, filtering inside the fractal due to characteristic func- 
tion depends only on the radius r, and this yields the 
integrations Q 

±,Sx{r)VrE,{r)dV = /^ x(/)[V.i?i(r')]r'^dr' 

= /;E.,ev„'5(r'-r,)G'(r')dr' (9) 

Here we define G{r) = r^Ei{r) and G'{r) = -^G. One 
obtains the integration of the electric field with the fractal 
density J^r-eVo ^^^' ~ ^j)- '^^^ latter corresponds to the 
fractal volume 

fj,(r) ^r" ^ r Yl '^(^' ~ r,y'^dr' . (10) 

•^0 r,eVo 



Following Ref. MT] (Theorem 3.1), we obtain 

^ G'ir')d^,ir') ^ -^J---J\r-r')^-'G'ir')dr' 



oir'G'ir). 



(11) 



Therefore, we consider the integration in Eq. © as 
the convolution integral with the averaged fractal den- 
sity (r - r')^-3. 

Now we estimate the integral /„ i?(r')Vr'x(r'')r' dr' . 
The fractal dust Vd at the Nth step of the construc- 
tion consists of balls B^ with the radius A^r. For 
example, A^v ^ Iq- In the limiting case one obtains 
Vd = limAT^oo U -^^ .1^]- ^^^ characteristic function 
for every ball is x(A7v) — 'd{r — rj) — Q{r — rj — Aiy). Dif- 
ferentiation of the characteristic function on the intervals 
[rj,rj+AN] yields Vrx(A7v) = S{r-rj)~S{r-rj-AN)- 



Therefore, for any interval An and at r 
with the electric field yields 



integration 



rll+A^E(ryWrX{r) - E{ry - E{r + AM){r + A^f . 

This expression is not zero in the limit A^v ^ 0. 
Let E{rj) is the electric field outside the ball -Bjv and 
E{rj + Ajv) denotes the internal electric field. The re- 
lation between them, due to Eq. (8.2) in Ref. lj| for 
polarization of a dielectric ball, is 



E{rj + An) == 



3£ri 



£,„(a;) +2ed 



E{r,) . 



Therefore, the shift for the electric field in the limit 
Atv -> is 



E{r,)-E{r,+Q)=E{r,)- 



n{L0) -Ed 

,{uj)+2ed 



(12) 



Finally, integration of the polarization charge term yields 

S;^E{r')Vr'X{r'y^dr' = oI^E{ryVrX{r) 

= IJ^T^^ Er.evo lo E{r')r"S{r' - r,)dr' . (13) 

Again, one obtains the integration of the electric field 
with the fractal density J2r ev ^(^' ^ ''i)- This corre- 
sponds to the fractal volume (ITUl) . and, hence, we con- 
sider the integration in Eq. (J13p as the convolution 
integral of Eq. (|TT|) with the averaged fractal density 
(r — r')^~^. Eventually, we obtain for the polarization 
charge term in Eq. ([7]) 



j;; E{r')\7r'X{r')r"dr' ^ 
^^j;;{r-rr-'Eir'y'dr' 
= o4^-2[G(r)+^or2]. 



(14) 



Here Eq = E • f is a projection of the external electric 
field on the radial direction inside the chosen spherical 
volume of the radius r and 



p{uj) 



Ed- emjuj) 
em(a;) + 2ed 



(15) 



Taking all these arguments and results in Eqs. pT|) 
and (fT4| . one presents Eq. Q in the coarse-graining 
form 



oir Gy) - qiioh^G'ir) - p(c.)o/,^-^G(r) 

2p{uj)Eo ^D ^g^ 



r(i:i + i) 



The Laplace transform can be applied to Eq. P^ . Since 
G{r = 0) = G"(r = 0) = (the electric field of frac- 
tal charge density diverges slowly than -^ @, [laDi one 
obtains for G{s) = C[G{r)] due to Eq. ©: 



G{s) 



2p{uj)Eo 



eds-' 



q{uj)s 



D-2 



-p{uj) 



(17) 



Note that the second term in Eq. p^ is q{Lj)G{r). 

Before arriving at the main result, let us consider the 
limiting cases. For £m(a;) = £d one obtains <?(w) = oo and 
p{io) = 0. This yields Ei{r) — 0, and the solution for the 
electric field is E{r) = Eq = Er- Another limiting case 
is |£m(w)| — )■ OO. In this case q{uj) — and p{uj) = —1, 
thus El ~ _Eo . Important result here is that permittivity 
of the mixture is approximately Ed that corresponds to 
the well known result in Ref. [IJ] [see Eq. (9.7)]. 



IV. SURFACE PLASMON RESONANCE 

The most interesting case is the SPR, when Re\e(jjj)\ = 
—2ed- This also known as the Frohlich resonance [l9ll20|. 
The permittivity of the fractal metallic nanostructure at 
the resonance condition is a complex value £m{'^) = £i + 
162, where £2/61 <^ 1 and is described by classical Drude 
formula 23] . At a small detuning from the resonance, 
when Re[e{Lu)] ~ ~2ed+S2 that corresponds to the width 
of the resonance in the frequency domain, p{uj) reaches 
the maximal values, that yields 

These expressions are inserted in Eq. p7)) . Before car- 
rying out the inverse Laplace transform £~^[G(s)], it is 
reasonable to simplify the second denominator. We re- 
cast Eq. (fT7|) in the form 



G{s) = 2p{uj)Eo J2 



[q{u)s^-^ + p{ujy 



qk+i 



(19) 



fc=0 



We take into account that for the scale r ^ 1, the Laplace 
parameter is small s ^ 1, and the binomial becomes 
approximately a monomial. For D < 2, the first term 
can be left, while for Z? > 2 the second term p{uj) is 
dominant. To keep the first term for D < 2, we have 
(1/s) ■^ r :^ {3ed/2e2)^-° ■ But this condition violates 



another restriction for r <C A/Zq, and p{uj) is the most 
important term in the binomial, and cannot be omitted. 
The Laplace inversion of Eq. (fT9|) can be performed using 
an expression for the Mittag-Lefher function [21 1 



£{^,0){zr'') 



27ri 



S-'-PgSr 



ds 



27r 



Z^ „^k+B ^* - Z^ 



k=0 



fc=0 



r{i^k + /3) 



(20) 



where C is a suitable contour of integration, starting and 
finishing at — oo and i/, /3 > 0. Comparing Eqs. (fT9| and 
([20|) . /3 = 4 and v = 1, one obtains for the electric field 



£;i(r) == 2p(w)i^or'£(i,4)b(w)r] 



(21) 



Since the argument of the Mittag-LefHcr function is large, 
its asymptotic behavior is 21] 



2EQr ,^y Eq 



ex — exp 






(22) 



Eventually, we arrived at the exponential (geometrical) 
enhancement and giant oscillations of the respond electric 
field due to the fractal geometry of the metal-dielectric 
composite. This result is valid for r > 1. At the same 
condition one can also choose the SPR at Reem(a;) = 
—2ed — TIS2, where ry is defined from the maximum en- 
hancement of the electric field. This yields q{uj) « 1/3 
and p{ui) w ed/s2- Thus the first term in the denomina- 
tor in Eq. (fTfl) can be neglected and the latter can be 
recast in the form 



Gis) 



2p{Lj)Eo 
q{u}) 



k=a 



,(fc+i)(r>-2)+3 



(23) 



This yields the following expression for the electric field 



Ei[r) = ^— £(D-2 

q{uj) 



3)hp(^)r^-7<z(^)]- (24) 



If 7] is chosen such that arg[ 



-p('^) 



< 



(_D-2)7r 



the 



L q(omega) -I ^ 2 

asymptotic behavior of the Mittag-Leffler function is ex- 
ponential '21] that leads to the giant geometrical en- 



hancement Ei{r) ^ exp 



(*)""■ 



that, obviously. 



takes place at the condition D > 2. 

When the resonance is exact and the detuning from the 
resonance is zero, then Re[p(a;)] = —1. The arguments of 
the Mittag-Leffler functions in both equations (PT|) and 
([M)) do not correspond to the exponential asymptotic 
behavior. In this case the geometrical enhancement of 
the electric field is absent, and the linear ("classical") 
enhancement [2Q| of the electric field takes place due to 
the imaginary part of the pre-factor: l'nip{uj) ~ ^ 3> 1. 



V. DISCUSSION ON THE GEOMETRICAL 
ENHANCEMENT OF THE ELECTRIC FIELD 

The obtained expressions in Eqs. (fT9|) and (|21|) (to- 
gether with developing the coarse grained Maxwell equa- 
tion in the form of Eqs. dlB]) and (HZ])) contain the main 
physical result. As seen, the resulting (enhanced) electric 
field depends on the parameter p{uj), which is the pre- 
factor and the argument of the Mittag-Leffler function 
(when r ^ 1) in Eqs. ([2T|) and (p4)) . and defined in Eq. 
([T5|) . This parameter is absorbtion efficiency in the elec- 
trostatic (quasi-static) approximation [20|. In our case, 
it describes polarization, and it is obtained in Eq. (J13p 
under evaluation of fractal boundary conditions. One 
should recognize that the linear (classical [20]) enhance- 
ment of the electric field is always takes place, due to Eq. 
([T^ that corresponds to the enhancement of the elec- 
tric field by a small particle and that is reflected by the 
pre-factor in Eq. ([T9|). For a fractal small composite of 
many particles the situation differs essentially. Here the 
geometrical enhancement of the electric field is due to 
the Mittag-Leffler function and it depends on the argu- 
ment of the complex value of p{ijj) . The exponential (ge- 
ometric) enhancement takes place when |p(w)] ^ 1 and 
Rep(w) > as in Eq. ([2T|), or Re[-p{uj)] > as in Eq. 
([Ml) . These conditions are fulfilled for those frequencies 
uj that are in the vicinity of the SPR: Rep(a;) = —2ed + A, 
where A ^^ £2 ^ £rf- 

We have to admit that the obtained expressions in the 
exponential forms are the upper bound of the electric 
field enhancement. The enhanced Ei does not exceed of 
the order of 10*(v/cm), otherwise the nonlinear effects 
become important that violates the linear quasi-static 
consideration. This restriction yields r{ed/2e2) < 20, 
which is a reasonable value for experimental realizations. 
Therefore, we have the light wavelength A ^ 10~*cm, 
fractal inhomogeneity size Iq ^ (10~^ -^ 10~^)cm, and 
£d/£2 '^ f^T '^ 5 -^ 15, where uj is the optical frequency, 
while T is the relaxation time. The latter value deter- 
mines Iq, as well, which was introduced above as a min- 
imum self-similarity size. Note that t < Tj = -2-, where 
Ts is the surface relaxation time, while vp is the velocity 
on the Fermi surface (for free electrons). Therefore, from 
the condition wr » 1, one obtains Iq ^ — ^ 10~^cm. 



Geometrical enhancement out of the SPR 



Equation (l^lT) describes geometrical enhancement of 
the electric field out of the SPR, as well. For those fre- 
quencies u! that p{uj) > 1, the argument of the Mittag- 
Leffler function is large, when 1 ^ r <C r" that leads 

(roughly) to the exponential Ei{r) ^ Eoexp[p{uj)r]. For 
example, when em(a;) = —Sd, the small imaginary part 
of em(w) is not important, and the condition for the 
large asymptotics of the Mittag-Leffler function is ful- 
filled. One obtains p{uj) w 2 and Ei (r) ^ Eq exp[2r] that 
yields an enhancement of the order of 10^ for r > 5. 



VI. CONCLUSION 

We developed an analytical approach, for description 
of the wave propagation in metal-dielectric nanostruc- 
tures in the quasi-static limit. The method is based on 
fractional calculus and permits to obtain an analytical 
expressions for the electric field enhancement. This ap- 
proach establishes a link between fractional geometry of 
the nanostructure and fractional integro-differentiation. 
An essential (geometrical) enhancement of the electric 
field is obtained for the surface plasmon resonance at a 
certain relation between permittivities of the host and 
fractal metallic nanostructure, when Re£m(a;) = —2ed 
with a suitable detuning. 

Important part of the analysis is developing convolu- 
tion integrals that makes it possible to treat the fractal 
structure. The initial Maxwell equation ^ is local, since 
Iq/X ^ 1 and space heterogeneity is accounted locally by 
virtue of the characteristic function. An accurate treat- 
ment of the fractal boundaries and recasting the Maxwell 
equation in the form of the convolution integrals by ac- 
counting fractal properties of the composite, eventually, 
leads to the coarse graining equations, which already take 
into account the space heterogeneity and nonlocal nature 
of the electric field and polarized dipole charges inside the 
composite. Therefore, the heterogeneity, caused by the 
fractal geometry, is reflected by the convolution of the 



averaged fractal density and the electric field, according 
Egs. pTj) and ([H]). It is necessary to admit that this 
transform from "local" quasi-electrostatics to the non- 
local, which takes into account space dispersion of per- 
mittivit y e (r) is mathematically justified and rigorous 
enough |17| . The obtained convolution integrals are aver- 
aged values, since the fractal density r^~^ is the averaged 
characteristics of the fractal structure, and, eventually, it 
determines the space dispersion of the permittivity e{r) 
of the mixture. 

Summarizing, we have to admit that observation of 
macroscopic Maxwell's equations is related to averaging 
of microscopic equations pJl . This procedure for fractal 
composite media is not well defined so far, since, accord- 
ing fractal's definition, averaging over an y fi nite volume 
depends on the size of this volume itself [I^ . The main 
idea to overcome this obstacle is to refuse the local prop- 
erties of equations and obtain nonlocal coarse graining 
Maxwell's equations, which are already averaged. 
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